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Abstract 

This paper studies asymptotic behavior of solutions of a free boundary problem modeling the 
growth of tumors with two species of cells: proliferating cells and quiecent cells. In previous literatures 
it has been proved that this problem has a unique stationary solution which is asymptotically stable 
in the limit case e = 0. In this paper we consider the more realistic case < e <C 1. In this case, 
after suitable reduction the model takes the form of a coupled system of a parabolic equation and 
a hyperbolic system, so that it is more difficult than the limit case e = 0. By using some unknown 
variable transform as well as the similarity transform technique developed in our previous work, we 
prove that the stationary solution is also asymptotically stable in the case < £ ^ 1. 
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1 Introduction 

It has long been observed that under a constant circumstance, a solid tumor will evolve into a 
dormant or stationary state [H 1271 128j . In a dormant state, the tumor does not change in size, while 
cells inside it are alive and keep undergoing the process of proliferation. In 1972 Greenspan established 
the first mathematical model in the form of a free boundary problem of a system of partial diff'erential 
equations to illustrate this phenomenon [26J . Since then an increasing number of tumor models in similar 
forms have been proposed by many different groups of researchers (see, e.g., [T], [5]-[S], [TB], [TS]-pD]. 
[3D]-[33], and the references cited therein). Rigorous mathematical analysis of such models has drawn 
great attention during the past twenty years, and many interesting results have been obtained, cf., [B] - 
[E], [m, EH - 123, [31], and references cited therein. 

In this paper we study the following free boundary problem modeling the growth of a solid tumor 
with two species of cells — proliferating cells and quiescent cells: 

ect = Ac~F{c) iorxenit), t > 0, (1.1) 

c = 1 for x e dn{t), t > 0, (1.2) 
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Pt + \/ ■{vp)^[KB{c)-KQ{c)]p + Kp{c)q for X e t > 0, (1.3) 

+ V • {vq) = Kq{c)p - [Kp{c) + KD{c)]q for x G t > 0, (1.4) 

p + g = l forxeO(i), t > 0, (1.5) 

^{vo\{n{t))) = [ V-v{x,t)dx for t>0. (1.6) 



Here Q{t) is the domain occupied by the tumor at time t, c = c(x,t), p = p{x,t) and q = q{x,t) are the 
concentration of nutrient, the density of prohferating cehs and the density of quiescent cells, respectively, 
and V = v{x, t) is the velocity of tumor cell movement. Besides, F{c) is the consumption rate of nutrient 
by tumor cells, Kb{c) is the birth rate of tumor cells, Kp{c) and Kq{c) are transferring rates of tumor 
cells from quiescent state to proliferating state and from proliferating state to quiescent state, respectively, 
and K]:){c) is the death rate of quiescent cells. 

Equation (1.1) describes diffusion of nutrient (regarded as one species) within the tumor, where 
£ = Tdiffusion/Tgrowth is the ratio of the nutrient diffusion time scale to the tumor growth (e.g., tumor 
doubling) time scale. Typically, Tdiffusion ~ Imin (see pp. 194-195 of [1 ), while Tg^owth ~ Iday, so 
that e <C 1. Equation (1.2) reflects the fact that the tumor receives a constant nutrient supply from 
its boundary, and we have rescaled c so that the constant supply of it from the tumor boundary is 
exactly 1. Equations (1.3) and (1.4) are conservation laws for the proliferating and quiescent tumor 
cells, respectively. Equation (1.5) reflects the fact that the mixture of proliferating and quiescent cells 
in the tumor has a constant density, and we have rescaled p and q such that this constant equals to 1. 
Finally, the equation (1.6) reflects the fact that the variance rate of the whole tumor volume is equal to 
accumulation of all local volume variance rate within the tumor (recall that V • v{x, t) is the local volume 
variance rate). We assume that F{c), Kb{c), Kp{c), Kq{c), Kpi{c) are C°° functions, and they satisfy 
the following conditions: 

(Ai) F'{c) > for c e R, and F(0) = 0. 

(A2) K'g{c) > and K'p{c) > for c e M, and ^^^(O) = Kp{0) ^ 0. 

(A3) K'^{c) < 0, K'q{c) < 0, Kd{c) > and Kq{c) > for c e R. 

(Ai) K'gic) + K'jj{c) > for c G R. 

(A5) K'p{c) + K'q{c) > for c e R. 

For simplicity we assume that the tumor is a strict spheroid of radius R{t), so that 

VL{t) ^ {x (^m? -.r ^\x\ < R{t)}. 

Accordingly, we assume that all unknown functions c, p and q are spherically symmetric in the space 
variable, i.e., 

c = c(r,t), p=p{r,t), q = q{r,t), 

and 

X 

V = v{r, t) — . 

r 

It follows that the above model reduces into the following system of equations: 

2 

ect(r, t) = Crr{r, t) + -Cr{r, t) - F{c{r, t)) for < r < R{t), t > 0, (1.7) 
r 

Cr{0,t)=0, c{R{t),t) = l fort>0, (1.8) 
Pt(r, t) + v{r, t)prir, t) ^f{c{r, t),p{r, t)) for < r < R{t), t > 0, (1.9) 
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where 



and 



2 

Vr{r,t) + -v{r,t) ^ g{c{r,t),p{r,t)) for < r < t > 0, (1.10) 

v{0,t) = for t > 0, (1.11) 
R{t) ^ v{R{t),t) for ^ > 0, (1.12) 

f{c,p)=Kp{c) + [KM{c)-KN(.c)]p-KMic)p^, 
9{c,p)=Km{c)p - Kd{c), 

Km {c)^ Kb{c)+Kd{c), Kn{c) = Kp{c) + Kgic). 



Global well-posedness of this model has been proved in the literature [TS] in a more general setting 
(see [13] for the limit case e = 0). In [T3] it is proved that under the assumptions {Aiy^A^) the above 
model has a unique stationary solution. In [SJ [5] it was proved that in the limit case e = 0, this unique 
stationary solution is asymptotically stable. The purpose of this paper is to prove that the result of 
[6j H] can be extended to the more realistic case £ ^ but small, namely, we shall prove that there exists 
eo > such that the stationary solution of the system (1.7)-(1.12) is also asymptotically stable in the 
case < e < eq. 

To present the precise statement of our main result, we make further reduction to the model as 
follows. Let 

c(f,i) =c(fe^(*\0, pif,t) ^p{fe''<-'\t), v{f,t) ^vife'^'\t)e^''-'\ i?(0 = e"(*\ (1.13) 

where < r < 1, t > 0. Then the system (1.7)-(1.12) is further reduced into the following system (for 
simplicity of the notation we remove all bars " ~ " ) : 

ee2^(*)ct(r,t) = c„(r,i) + +ee2^(*Vu(l,t)]c^(r,i) - e2'=(*)i^(c(r,<)) for < r < 1, t > 0, (1.14) 

Cr{0,t)=0, c(l,i) = l fort>0, (1.15) 

Pt{r,t)+w{r,t)pr{r,t)^f{c{r,t),p{r,t)) for < r < 1, i > 0, (1.16) 
2 

Vr{r,t) + -v(r,t) = g(c(r,t),p(r,t)) for < r < 1, t > 0, (1.17) 
r 

v{0,t) = for t > 0, (1.18) 

z{t) = w(l,i) for t > 0, (1.19) 

w{r, t) = v{r, t) ~ to(1, t). (1.20) 



where 



Let (c*(r),p*(r),w*(r),z*) be the stationary solution of the system (1.14)-(1.20) ensured by [13]. 
The main result of this paper is the following: 

Theorem 1.1 Let the assumptions {Ai)-{A^) he satisfied. There exist positive constants £q and pL* 
such that for any < £ < eq and < /i < /i*, there exist corresponding positive constants 6 and C such 
that for any time- dependent solution {c{r,t),p{r,t),v{r,t), z{t)) of the system (1.14)-(1.20), if the initial 
data {co{r),pQ{r),zo) = (c(r, 0),p(r, 0), z(0)) satisfy 

co{r), c'o{r)eC[0,ll €^,(0) = 0, co(l) = 1, < co(r) < 1 /or < r < 1, (1.21) 

po{r), r{l-r)p'o{r) eC[OA], po(l) = 1, < po(r-) < 1 forO<r<l, (1.22) 
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max \co{r) - c*(r)| < 6, sup |co(r) - c* (r)| < 5, (1.23) 

0<'-<l 0<r<l 

ma.x \po{r) - p*{r)\ < 6, snp r{l-r)\p'g{r) - p*' {r)\ < 6, and \zo-z*\<S, (1.24) 

0<r<l 0<T<1 

then for all t > we have 

max [\c{r,t)-c*(r)\+ max \cr{r,t) - c*' {r)\ + max \ct(r,t)\] < CSe'''^ (1.25) 

0<r<l 0<r<l 0<r<l 

max \p{r,t) -p*ir)\ + max r{l-r)\pr{r,t) - p*' {r)\ + max \pt{r,t)\] < C6e'~'^\ (1-26) 

0<r<l 0<r<l 0<r<l 

and 

\zit)-z4 + \z{t)\<C6e~''\ (1.27) 

Remark: Readers who are familiar with the hteratures [23], [T^] and [7] might wish to use the 
iteration technique used there to prove the above theorem. The author is not successful in trying to do 
so. The obstacle lies in making estimates for |f — S{f, t, s)\ and related quantities, where S{f, t, s) denotes 
the similarity transform introduced in [5] (see (4.16) in Section 4 below). In this paper we use a different 
approach; see Section 2 for illustration. 

The structure of the rest part is as follows: In the next section we make further reduction to the 
system (1.14)-(1.20) and illustrate the main idea of the proof of the above theorem. In Sections 3-5 we 
make preparations for the proof of Theorem 1.1, which is given in the last section. 



2 Reduction of the problem 



We define a function m in [0, 1] x M as follows: For any z G M, the function c(r) — m{r; z) {0 < r < 1) 
is the unique solution of the following boundary value problem: 

2 



c"(r) + -c'(r) = e^^F{c{r)) for < r < 1, 



r 



(2.1) 



c'(0) = 0, c(l) = l. 



SinceFe C°° (M) , we have m e C°°([0, 1] xR). For the solution (c(r, t),p(r, t), ^(r, i), z(t)) of the problem 
(1.14)-(1.19), we let 

ri{r,t) = c{r,t) - m{r;z{t)). (2.2) 
It can be easily seen that rj satisfies the following equations: 

ee2^(*)77t(r, t) = 77,,(r, t) + + ee^'^'^rv{l, t)] 77,(r, t) - e^'^*'>a{r; Tj{r, t), z{t))Tj{r, t) 

-se^^i^)v{\,t)[m,Xr;z{t)) -rmr{r]z{t))\ for < r < 1, i > 0, (2.3) 
7^(0, t) = 0, ?7(l,t) = for i > 0, 
where a is a function in [0, 1] x M x M defined as follows: For r e [0, 1] and y, z e R, 

a{r]y,z)= [ F'{m{r- z) + ey)de. (2.4) 
Jo 

Conversely, if r]{r,t) satisfies (2.3) and z{t) satisfies (1.19), then c{r,t) = m{r; z{t)) + ri(r,t) satisfies (1.14) 
and (1.15). Hence, under the transformation of unknown variables 

{c{r,t),p{r,t),v{r,t),z{t)) ^ {i^{r,t),p{r,t),v{r,t), z{t)) (2.5) 
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given by (2.2), the system (1.14)-(1.20) is equivalent to the system (1.16)-(1.20) (with c{r,t) replaced by 
m(r; z{t)) + r]{r, t)) coupled with (2.3). 

We introduce three maps v,w,f : C[0, 1] x C[0, 1] x M — >■ C[0, 1] and a functional g in C[0, 1] x 
C[0, 1] X R as follows: For r],pe C[0, 1] and zeR, 



w{r;r],p,z) 



9{'m{p,z) + r]{p),p{p))p'^dp for < r < 1, 

0, for r = 0; 
w{r;r],p,z) =w{r;r],p,z) - rv{l;r],p, z), 

{{r;T],p, z) =f{m{r,z) + r]{r),p{r)), 

&{V,P, z) =v(l; r),p, z)= I g{m{r, z) + ri{r),p{r))r'^dr. 

Jo 

It can be easily seen that with t] as in (2.2), we have 

v{r, i)=v(r; r?(-, t),p{-,t), z{t)), 
w{r, t)=w{r; r]{-,t),p{-,t), z{t)), 
f{c{r, t),p{r, t))=f (r; rj{;t),p{;t), z{t)), 
v{l,t)=M,t),p{;t),z{t)). 

Hence, under the transformation of unknown variables given by (2.2) and (2.5), the system (1.14)-(1.20) 
reduces into the following system: 



ee^'=(*)ryt(r,t) = r]rrir,t) + [- + ee'^^''*^rg{r],p, z)]r]r{r,t) - e^'-^'^'>a{r]rj,z)r] 
—ee'^^^*^g{ri,p, z)[m,z{r; z) ~ rmr{r; z)] for < r < 1, t > 0, 
r]r{0,t) = 0, r]{l,t) =0 for t > 0, 
Pt + w{r;ri,p,z)pr ^ f{r;r],p,z) for < r < 1, t > 0, 
z = g(ry,p, z) for t > 0. 

We shall treat the above system in the following way: Let 

y = jr; = r]{r, t) G C([0, 1] x [0, oo)) : r]{r, t) is differentiable in r, r]r G C([0, 1] x [0, cx))), 

r?r(0,t) = 7?(l,t) = 0, > 0, and ||r?||r sup e''*|r?(r,t)| + sup e*** r/r(r-,i) < ooj, 

0<T-<1 0<r<l J 



(2.6) 



0<T- 

t>0 



0<r 

i>0 



where /i is a positive constant to be specified later. It is clear that (Y, || • \\y) is a Banach space. Let 5 
and 5' be two positive constants to be specified later. We assume that the initial value {co{r) , po{r) , zq) 
of {c{r,t),p{r,t), z{t)) satisfies the conditions in Theorem 1.1. Given r] gY satisfying the conditions 



< S' and ri{r, 0) = r]o{r) =^ co(r) - m(r, zq), 
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we first solve the initial value problem 

f 

pt + w(r; r],p, z)pr = f (r; r],p, z) for < r < 1 and t>Q, 

z~g{ri,p,z) for t > 0, 

^p{r, 0) = po{r) for < r < 1, and ^(0) = zo, 

and next solve the following initial-boundary value problem: 
' 2 

ee^^Wrjt = r?rr + [- + ee'^''^*\s{r],p, z)]fir - e^''^''^a{r;r], z)fj 

—£e^^(.t)g^(^'q^p^z)[mz{r;z) — rmr{r;z)\ for < r < 1, t > 0, 

fir{0,t)=0, fj{l,t) = fovt>0, 

fi{r, 0) = r]Q{r) for < r < 1. 



(2.7) 



(2.8) 



We then obtain a map r] ^ fj from the closed ball B^i (0) in Y into Y. We shall prove that there exists 

£0 > such that for any Q < e < £i^/ii 5 and 5' are sufficiently small then this is a self-mapping in ^^'(O). 
Moreover, we shall use the Schauder fixed point theorem to prove it has a fixed point in this ball. The 
desired assertion then follows. 



3 A preliminary lemma 

Lemma 3.1 Assume that F'{c) > for c e R and F{0) = 0. For all < r < I and z gM. we have: 
(1) < m(r; z) < 1; (2) mr{r; z) > 0; (3) m^(r; z) < 0; 



(4) mrir;z)<^rF{l)e^'; 



(5) m.{r;z)>-^{l-r^)Fil)e'\ 



(6) --F(l)e^^ < mrr{r;z) < F{l)e''; 



(7) m,{r;z)>--{l-r^)F{l)e^\ 

Proof. The assertions (1)~(3) are well-known. Let c(r) = m{r;z). From the first equation in (2.1) 
we have 

(c'(r)r2)' = e^'-F{c{r)y. 
Integrating both sides from to r and using the assertion (1) we obtain 



c'(r) 



o2z 



J^^ F{c{p))p'dp<^rF{l)e'^. 



(3.1) 



This proves the assertion (4). Next, differentiating all three equations in (2.1) in z we see that the 
function Cz{r) = mz{r; z) is the solution of the following boundary value problem: 

c''(r) + -c' (r) = e^' F' {c{r))cz{r) + 2e^''F{c{r)) for < r < 1, 
r 

^4(0) =0, c,(l) = 0. 
Since F'{c) > 0, Cz{r) < and c(r) < 1, we see that 

(4(r)r2)' = e^'F'{c{r))cz{ry +2e^''F{c{r)y < 2e=^^F(l)r^ 
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It follows that 

4(r) < ^F(l) rp^dp < lrF{l)e'\ 
'''Jo -J 

Hence 

-c.(r) < lF{l)e'^ J' pdp = 1(1 - r2)F(l)e2^ 
This proves the assertion (5) and completes the proof of Lemma 3.1. □ 

4 Decay estimates of the solution of (2.7) 

We denote by Cy[0, 1] the following function space: 

Ci[0, 1] = e C[0, 1] n C\0, 1) : r(l - r)</)'(r) e C[0, 1]}, 

with norm 

||(/)||ci[o,i] = max |0(r)| + sup |r(l - r)<j>'{r)\ for e Ci[0, 1]. 

0<r<l 0<r<l 

It is clear that this is a Banach space. Next we denote 

X = C[0, 1] X K and Xq = C^[0, 1] x K; 

they arc both Banach spaces. Given t] G Y, we introduce a map F : [0, oo) x Xo ^ X as follows: For 

f > and [/ = {p,z) G Xq, 

¥(t,U) = (-w(-;r/(-,t),p,z)p'+f(-;?7(-,t),p,z),g(r/(-,t),p,z)). 

Then the system (2.7) can be rewritten as the following differential equation in the Banach space X 
(regarding Xq as a subspace of X and F(-, as an unbounded nonlinear operator in X with domain Xq): 

U = ¥{t, U) for t > 0. (4.1) 

It is clear that [t ^ F(<, •)] € C([0, oo), C°°(Xo, X)), and 

F(t,C/) = Ao(t,C/)J7 + Fo(t,?7) for [/ e Xq, t > 0, 

where 

Ao(t,C/)V'=(-w(-;r?(-,t),p,^)g',0) for (7 = (p, z) e X, F = (g,y)eXo, i > 0, 
Fo(t,C/) =(f (•;//(•, t),p,z),g(r?(-,i),P,^)) for C/ = (p, z) G X, t > 0. 

It is also clear that [t ^ Ao{t, •)] S C([0,(X)),C°°(X,.if(Xo,X))) and [t ^ ¥{t, •)] e C([0,(X)), 
C°°(X,X)). Let [/* = ip*,z*) andV = U- U*. Then 

F(i,?7) =Ao{t,U)U + ¥o{t, U) 

=Ao{t, U* + V){U* + V)+ ¥o{t, U* + V) 

=Ao(t, u* + y)y + Ao(t, [/* + F)c/* + ¥o{t, u* + y) 

={Ao(t, [/* + y)t/ + [duAl{U*)V]U* + 9r/Fa(;7*)F} 
+{[at/Ao(i, [/*)F - duAl{U*)V]U* + [%Fo(t, C/*)F - du¥l(U*)V]] 
+{[Ao(t, [/* + y) - Ao(t, U*) - duAoit, U*)V]U* + [Fo(t, U* + V) 
-¥o{t, U*) - du¥o{t, U*)V]} + {Ao{t, U*)U* + ¥o{t, U*)} 

=Ait, V)V + M{t)V + Git, V) + G{t), 
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where duAo{t,U), du¥Q{t,U), duA*{U), du¥*{U) denote the Frechet derivatives of Ao{t,U), ¥o{t,U), 
Aq{U), Vq{U) in the variable U, respectively, where 



*o{U)V ={-w{-,0,p,z)q',0) ^ \im Ao{t,U)V ior U = {p, z) e X, V ^ {q,y) e Xo 
F5(C/) =(£(•; 0,p,z),g(0,p, z)) = lini Fo(i,[/) for U = {p,z)eX, 



t^OC 

In 

t^OO 

and ioTVeX,We Xq and t > 0, 



A{t, V)W =Ao(i, U* + V)W + [duAo{t, U*)W]U* + du¥o{t, U*)W, 

M{t)V ^[duAoit, U*)V - duA*iU*)V]U* + [duFoit, U*)V - du¥*{U*)V], 
<G{t, V) =[Ao(i, U* +V)- Ao{t, U*) - duAo{t, U*)V]U* + [¥o{t, U* + V) 
-¥o{t,U*) ^ du¥oit,U*)Vl 
G{t) ^Ao{t, U*)U* + Fo(t, U*) = F(t, [/*). 

Hence (4.1) can be rewritten as follows: 

V ^ A{t,V)V + M{t)V + <&{t,V) + G{t) ior t>Q. (4.2) 

Note that the above equation is quasi-linear, and 

A{t,V) = Ao(t,C/* + F)+B*, 

where 

M*W = [duAl{U*)W]U* + du¥l{U*)W for W & X. 

It is easy to see that B* e .^{X) n ^{Xq) (cf. Corollary 3.2 of [8 ). Moreover, using the mean value 
theorem we can easily prove that B e C([0, oo), ^(X) n ^{Xq)), and there exists constant C > 
independent of rj such that 

WmW^iX) + mm^ixo) < Clmax |r;(r,t)| + max \vr{rM < CSe-^^* (4.3) 

0<r<l 0<r<l 

for all t > 0. Besides, since rj and rjr are bounded functions. It is also easy to see that there is a positive 
constant C independent of t] such that 

\m,V)\\x<C\\V\\l, (4.4) 

\m,V)\\x„ < C\\V\\l^, (4.5) 
||G(i,Fi) - G(t, V2)\\x < C{\\Vi\\x + \\V2\\x)\\Vi - V2\\x (4.6) 
mt,V^)~G{t,V2)\\xo<Cm\\x„ + \\V2\\xo)\\Vi-V2\\xo (4.7) 

for all i > and small HV^Hx, ||Vi||xoj ll^2||xo (cf- Corollary 3.3 of [8]). As for the last term on the 
right-hand side of (4.2), we have: 

Lemma 4.1 Let rj Cz Y and \\'I]\\y < S. Then G G C([0, oo), X^) and there exists a positive constant 
G indepedent of 5 such that 

\\G{t)\\x^ < G5e~''^ for t > 0. (4.8) 

Proof: The assertion that G G G{[0,oo), Xq) follows immediately from Lemma 3.1 of j8j. To prove 
(4.8) we note that 

G(t) = F(t, U*)={~ w(-; r/(-, t),p\z*)p*' + f{--T^{-,t),p\z*), g{ri{;t),p*,z*)). 
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Since (c*,p*, u*, z*) is stationary solution of (1.14)-(1.20), we have 

c*{r) = m{r,z*), -~v*{r)p*' {r) + f{c*{r),p*{r)) ^ 0, and 

v*{l)= f\{c*{s),p*{s))s^ds^O. 
Jo 

Moreover, since v*{l) = 0, we have 

\w{r;r^{-,t),p*,z*)-v*{r)\ <\^{r;7j{;t),p* , z*) - v*{r)\ + r\vil;v{;t),p* , z*) ~ v*{l)\ 

<\ f \g{c*{p)+q{p,t)y{p))-g{c*{p),p*{p))\p^dp 

Jo 

+r [ \gic*{p)+,^{p,t),p*{p))-g{c*ip),p*ip))\p^dp 
Jo 

<Cr max \ri{r,t)\ for < r < 1, t>0. 

0<r<l 

Hence 

||G(t)|lx=max \wir;r,{;t),p* , z*)p*' (r) ^ f{-,T^{;t),p* , z*)\ + \g{ri{;t),p* , z*)\ 

0<r<l 

< max r^^\w(r; ri(- t), p* , z*) — v*(r)\ ■ max rp* (r) 

0<r<l 0<r<l 

+ max \f{c*{r)+rj{r,t),p*{r)) ~ f {c* (r) , p* {r))\ + \v{l;f^{-,t),p*,z*)~v*{l)\ 

0<r<l 

<C max \Tj{r,t)\ < CSe^^'^ for t > 0. 

0<r<l 

Similarly, by using the above-mentioned equations for {c* ,p* , u*, z*) and the fact that rp* (r), 
r^p* (r) e C[0, 1] (cf. Lemma 3.1 of [S]) we can also prove that 

max |r(l - r)|-[w(r; T^i;t),p* , z*)p*' (r) - f (•; 77(-, t),p*,z*)]\ 
0<r<l or 

<C max \ri(r,t)\ +C max \r]r(r,t)\ < CJe"''* for t > 0. 

~ 0<r<l 0<r<l I ' \ ' 'I — — 

Hence (4.8) holds. □ 

Let V S C([0, oo),X) be given. By using some similar arguments as in the proof of Lemma 4.2 of 
[5^, we can easily prove that the family of unbounded linear operators {A(i, V) : i > 0} in X is a stable 
family of infinitesimal generators of Cq semigroups in X, and {A{t, V) : t > 0}, their parts in Xq, is a 
stable family of infinitesimal generators of Co semigroups in Xq- By Theorem 3.1 in Chapter 5 of [29] . 
it follows that A(i, V) {t > 0) generates an evolution system U(t, s\V) (t>s>Q). Moreover, by using 
a similar argument as in the proof of Lemma 4.3 of [5], we can easily prove that for any i7o G and 
s > 0, the initial value problem 

U{t) = A{t,V)U{t) (for t > s) and U{s) = Uq (4.9) 

has a unique solution U = U{t; s, Uq) £ C([0, oo), Xq) f| C^([0, oo), X). By Theorem 4.2 in Chapter 5 of 
[29] . it follows that U{t; s, Uq) = l]{t, s; V)Uo for aU t > s > 0, and, furthermore, for any F e C([0, oo), Xq) 
the expression 

U{t) = U(i, 0; V)Uo + I l]{t, s: V)F{s)ds for t > 0. 

"'0 
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gives a unique solution U G C{[0,oo),Xo)f]C^{[0,oo),X) of the initial value problem 

U{t) = A{t, V)U{t) + F{t) (for t > 0) and U{0) = Uq. 

Hence, given Vq G Xq, the unique solution of the equation (4.2) subject to the initial condition 1^(0) = Vq 
is also a solution of the integral equation 

V{t)=U{t,0;V)Vo+ [ i]{t,s;V)[M{s)V + G{s,V)ds + G{s)]ds (4.10) 

and vice versa, i.e., the equation (4.2) subject to the initial condition V{0) = Vq is equivalent to the 

above integral equation. 

Lemma 4.2 Let rj G C([0, 1] x [0,oo)) and V = V{t) G C([0,oo),X) be given such that r]r G 
C([0,1] X [0,oo)), 

sup \v{r,t)\ + sup \Vr{r,t)\ <Se-^'* and \\Vit)\\x < CoSe-^"* fort>0, (4.11) 

0<r<l 0<r<l 

where Co, S and fi are given positive constants. There exists a positive constant fi* determined by the 
operator A*(0) — lim A(t,U*) (therefore it is independent of the constants Co, S and fi), such that for 

any < /x' < /x* there exists corresponding constant i5o > (depending on Cq, h and fi') such that if 
< ^ < ^0 then the following estimates hold: 

||U(t, s, V)\\l(x) < Cie-'^' /or i > s > 0, (4.12) 

mt, s, V)\\Lix,) < Cae-'^'^*-^) fort>s>Q, (4.13) 
where C\ and C2 are positive constants depending only on fj,' and not depending on 11 and Co . 

Proof Let V(t) = (ip(-,t),C(t)) and set 

p(r, t) = p* (r) + ip(r, t) , z(t) = z* + C(t) . 

We denote 

w(r,t) =w(r;r](-,t),p(-,t),z(t)) 

=K [ 9(m{p,z(t)) + 'q(p,t),p(p,t))p'^dp-r I g(m(p, z(t)) + 'q(p,t),p(p,t))p'^dp 
7" Jo Jo 

for < r < 1, t > and w(0,t) = for t > 0. Similarly as in the proof of Lemma 4.1 we can prove the 
following estimates: 

max \w(r,t) - v*(r)\ < C5r(l - r)e-''* for t > 0, (4.14) 

0<r<l 

dw(r,t) ^*'^^-) <^jg-Mt fort>0. (4.15) 



max 

0<r<l 



dr 



Indeed, for < r < 1/2 we have 

\w(r,t)-v*(r)\ <\ [ \g(m(p,z(t))+'n(p,t),p(p,t))-g(m(p,z*),p*(p))\p''dp 
^ Jo 

+r [ \g(m(p,z(t))+r](p,t),p(p,t)) - g(m(p,z*),p*(p))\p'^dp 
Jo 



<Cr[\z(t) - z*\ + max \7^(r,t)\ + max \p(r,t) - p*(r)\] 

0<r<l 0<r<l 

<CrSe-i'^ for t > 0, 
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and for 1/2 < r < 1 we have 

\w{r,t)-v*{r)\ <{^-r) j \g{m{p, z{t)) + rj{p,t),p(p,t)) - g{m{p, z*),p* {p))\p^dp 



\g{m{p, z{t)) + ij{p,t),p{p,t)) - g{m{p,z*),p*{p))\p^dp 



<C{1 - r)[\z{t) - 2*1 + max |77(r, t)| + max \p{r,t) - p* {r)\] 

0<r<l 0<r<l 

<C(1 - r)Se-^'* for t > 0. 

Combining the above two estimates together we obtain (4.14). The proof of (4.15) is simplar; we omit it. 

Having proved (4.14) and (4.15), it follows that all the results in Section 5 of [8 can be applied to the 
family of similarity transforms f t-^ r = S{f; t, s) {t > s > 0) of the unit interval [0, 1] to itself obtained 
from solving the initial value problem 

' dr dr 

+ v*ir)i^ ^ w{r,t) for < r < 1, t > s, 
dt dr (416) 

r\t=s — r for < f < 1- 

The estimates (4.12) and (4.13) then follows from a similar argument as in the proof of Lemma 6.4 of [8]. 
This proves Lemma 4.2. □ 

Let p* be the constant as in Lemma 4.2 and fix a number < p < p*. Later on the notation p 
always denotes this fixed number. 

Lemma 4.3 Let rj E C([0,1] x [0,oo)) be given such that rjr E C([0,1] x [0,oo)) and 

sup \r){r,t)\ + sup |?7r(r, t)| < (Je"'"* fort>0, (4.17) 

0<r<l 0<r<l 

where d and p are given positive constants, S sufficiently small and p < p* . Let Vq G Xq be such that 
ll^llxo ^ ^- Under these conditions, the equation (4.2) subject to the initial condition V{0) = Vq has a 
unique solution V E C{[0,oo),Xq) f] C^{[0,oo), X) satisfying the following estimates: 

\\V{t)\\x„ < CSe-^\ \\V{t)\\x < C5e-^* for t > 0, (4.18) 

where C is a positive constant independent of rj, Vq and 5. 

Proof: Let M be the set of all functions V = V{t) E C([0, co),X) satisfying the following conditions: 

V^(0) = yo, \\V {t)\\x <C5e->'^ loi t>0, (4.19) 

where C is a positive constant to be specified later. We introduce a metric d on M by defining d{Vi, V2) = 
supj>Q e***!! Vi(i) — V2(i)||x for ¥±,¥2 E M. It is clear that (M, d) is a complete metric space. In what 
follows we split into two steps to prove that the equation (4.2) subject to the initial condition 1^(0) = Vq 
has a unique solution in M provided S is sufficiently small, and the solution also satisfies the other two 
estimates in (4.19). 

Step 1: We prove that if 5 is sufficiently small then for any V E'M., the initial value problem 

U{t) = k{t, V{t))U{t) + M{t)U{t) + G{U{t)) + G{t) for t > 0, 

(4.20) 

U{0) = Fo- 



il 



has a unique solution U e C([0,oo),Xo) fl C^([0,oo),X) satisfying the following estimates: 

\\U{t)\\x<C5e-^'\ \\U{t)\\x,<C'5e-^'\ \\U{t)\\x < C de-^'' for t > 0, (4.21) 

where C is the same constant as that appears in (4.19), and C is another positive constant to be specified 
later. To this end we let 

M = {[/ e C([0,oo),Xo) : \\U{t)\\x < Cbe'^^ and ||C/(t)|Uo < Cbe'^^ for t > 0}, 

and introduce a metric d on it by defining d{Ui,U2) = supj>Q e''*||i7i(i) — U2{t)\\xo for Ui,U2 S M. 
(M, d) is clearly a complete metric space. Given U G M, we consider the following initial value problem: 

= A{t, V{t))U{t) + M{t)U{t) + G{U{t)) + G{t) for t > 0, 

Since U{t) e C([0, oo), Xq), we have <G{U{t)) e C([0, oo), Xq). Since also G € C([0, oo), Xq), it follows 
that the above problem has a unique solution U e C([0, oo), Xq) fl C^([0, oo), X), given by 



U{t) 



\i{t,Q,V)Vo+ [ l]{t,s,V)[M{s)U{s) + G{U{s))ds + G{s)]ds. (4.22) 
Jo 



Choose a number fj,' such that fi < fi' < ii* . By Lemma 4.2, if 5 is sufficiently small then there exist 
constants Ci,C2 > depending only on /x' such that for any V = V{t) e C([0,oo),X) satisfying (4.19), 
the following estimates hold: 

l|U(t, s, V)\\l(x) < Cie-'''^*-'*) for t > s > 0, (4.23) 

||U(t, s, y)llL(Xo) < C2e-'''(*-«) for t > s > 0, (4.24) 

Using (4.22), (4.23), (4.4), (4.8), the fact that \\Vo\\x < \\Vo\\xo < S and the condition ||?7(i)||x < C6e-^'* 
we see that for some positive constant C", 



\\U{t)\\x<C,e-''''\\Vo\\x +C, f e-'''(*-«)[||B(s)[/(s)||x + \mU{s))\\x + \\G{s)\\x]ds 

t 

<Ci6e-^'' + CiC" f e-^'(*-^) [CS^e-^"' + C^S^e-^^' + Se'^^jds 
° t „ t 

<(l + 2C"C^6)c,S [ e-'^^'-'^e-^^'ds + ^^f-^e-'^'' I e^^'-^^'ds 
^ ^ \ M-M Jo 

<(l + -C"CH+-^\c^5e-''K 

\ jJL jJL' — llJ 

( C" \ 

Hence, if we first choose the constant C > sufficiently large such that C > 2 H Ci and next 

V /i' — /X/ 

2 

choose (5o > sufficiently small such that —C'C^do < 1, then for all < J < (5o we have 

M 

\\U{t)\\x < CSe-^"* for all i > 0. 

Similalrly, by using (4.22), (4.24), (4.5), (4.8), the fact that \\Vo\\xo < S and the condition ||?7(i)||xo < 
C"(5e~^* wc sec that by first choosing C sufficiently large and next choosing further small (when 
necessary), we also have 

\\Uit)\\xo < C'5e->'^ for all t > 0. 
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Hence U G M. We now define a map S : M — > M by setting S(U) = U for every U G M. A similar 
argument by using (4.7) instead of either (4.4) or (4.5) we can prove that by choosing So further small 
when necessary, S is a contraction mapping. It follows from the Banach fixed point theorem that S has 
a fixed point in M, which is clearly a solution of the problem (4.20) in C([0, oo), Xq). Uniqueness of the 
solution follows from a standard argument. The assertion that the solution U G C^{[0,oo), Xq) and it 
satisfies the third estimate in (4.21) are easy consequences of the equation in the first line of (4.20) (cf. 
the proof of Lemma 7.1 of [5]). 

The assertion obtained in Step 1 in particular implies that for every V in'M., the solution U of (7.2) 
also belongs to M. Thus we can define a mapping S : M — > M as follows: For any F S M, 

S{V) = U = the solution of (4.20). 

Step 2: We prove that if S is sufficiently small, S is a contraction mapping. For this purpose, let 
Vi,V2eM and denote J7i = S(Vi), C/a = 8(^2) and W = Ui~U2. Then W satisfies 

( dW^ ^ V,{t))W{t) + [Ait, Fi(t)) - A(t, V2mU2{t) 
dt 

+M{t)W{t) + [G{Ui{t)) - G{U2{t))] for t > 0, 
[VF(0) = 0, 
so that 

W{t) ^ f l]{t,s,Vi){[A{s,Vi{s)) - A(s,V2{s))]U2{s)ds + M{s)W{s) + [<G{Ui{.s)) - G(U2{s))]}ds. 
Jo 

It can be easily shown that (cf. the proof of Lemma 7.2 of [8]) 

\ms,Vi{s)) ~ A{s,V2{s))]U2{s)\\x < C\\Vi{s) ~ V2{s)\\x\\U2{s)\\xo < CSe-^f^'di^^ (4.25) 
Besides, from (4.3) we have 

\ms)W{s)\\x < C6e-'''\\W{s)\\x < C6e-^^-'d{Ui,U2), 

and from (4.6) we have 

||G([/i(s)) - G{U2{s))\\x<C{\\Ui{s)\\x + \\U2{s)\\x)\\Ui{s) - U2{s)\\x 
<C(5e-''^||T^(s)||x < CSe-^'"'d{Ui,U2). 

From these relations and (4.21) we get 

\\Ui{t) -U2{t)\\x<CSd{Vi,V2) f e^''(*~'*)e-2f-ds + CM(C/i,C/2) / er^''^'-'^ e'^^'^'ds 

Jo Jo 

<CSe-^'*d{Vl, V2) + C(5e-^*d([/i, f/2), 

which yields d{Ui, U2) < C^d(Vi, V2) + CSd{Ui, 1/2). The desired assertion now easily follows. 

It follows that if 6 is sufficiently small then the map S has a fixed point ?7 in M. Since the image 
of S is contained in M, we obtain the assertion of Lemma 4.3. This completes the proof of Lemma 4.3. 
□ 

Remark: A similar argument as in Step 2 of the above proof shows that if Vi , V2 are solutions of the 
equation (4.2) with respect to 771, 772 € K respectively (with same initial data), then we have the following 
estimate: 

d{VuV2)<CS sup e'''\mir,t) - 7j2ir,t)\. (4.26) 

0<r<l 
f >0 
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Indeed, to emphasize the dependence of the operator A{t,V) on rj we redenote it as A{r],V). Similarly 
as in (4.25) we have 

||[A(r?i(., s), V^{s)) - A{rj2{; s), V2is))]U2{s)\\x 
<C sup \mir,s)-m{r,s)\\\U2is)\\xo+C\\V^{s)-V2{s)\\x\\U2is)\\xo 

0<r<l 

<Ce-^^'6 sup e'''\riiir,t)-7j2{r,t)\+C6e-^'''d{Vi,V2). 

0<r<l 
t>0 

Prom this inequality and a similar argument as before we obtain the desired assertion. 
The estimate (4.26) will be useful in Section 6. 

5 Decay estimates of the solution of (2.8) 

Lemma 5.1 Let {r],p,z) = {r]{r,t),p{r,t), z{t)) (0 < r < l,t > 0) be given such that r],r]r G 
C([0,1] X [0,oo))(= C([0,oo),X)), p e C([0,oo),Xo), z e Ci[0,oo). Assume that 

sup \r]{r,t)\ + sup \r]r{r,t)\ < de'i^* for t > 0, (5.1) 

0<r<l 0<r<l 

sup \p{r,t) -p*{r)\ + sup |r(l - r)\pr{r,t) - p*' {r)]\ < ^e"'^* for t > 0, (5.2) 

0<r<l 0<r<l 

\z{t)- z*\ + \z{t)\<6e-^'* fort>0. (5.3) 
Let fj = fi{r,t) {0 < r < l,t > 0) be the solution of (2.8) with initial data fi{r,0) = r]o{r) satisfying 

sup |7?o(r-)| + sup |r7o(r)| < s'6, (5.4) 

0<r<l 0<r<l 

where e' is a given small positive constant. Under these conditions, there exist constants eo,£o > 
C > independent of 6 such that if < s < eq and < s' < s'q then 

\fl{r,t)\ + \fir{r,t)\ < C{e + e')6e-''' (5.5) 

for alio <r <1 and t>0. 

Proof Let co = min F'{c) > and A = coe2(^' Since s{0,p*,z*) = v*{l) = 0, we have 

— 5o<c<l+i5o 

Siv{;t),p{;t),z{t))=\s{v{;t),p{;t),z{t)) - s{0,p*,z*)\ 

< I \9{m{p,z{t)) + r]{p,t),p{p,t)) - g{m{p,z*),p*{p))\p^dp 
Jo 

<C[ max \r]{r,t)\ + max |p(r,t) - p*{r)\ + |^(i) - z*\] 

0<r<l 0<r<l 



Besides, we also have 



<CSe-i^\ (5.7) 

z* -6o<z* - (5e-^* < z{t) <z* + Je"^* < z* + 6o, 
\m;,{r; z) - rnirir; z)\ < Ce^^*) < Ce'^''^'+^°\ 



a{r; r]{r, t), z{t)) = / F'(m(r, z{t)) + er]{r, t))de > cq. 
Jo 
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It follows that 

|£e2^(*)g(j?,p, z) [m,{r; z) - rmr{r; z)] \ < Ce^e"^*, 

and 

e2^Wa(r;r?,z) > coe^^''-^"^ = X. 
Using these estimates, we can easily verify that 



{r,z) = 2X-^CeSe-t'' + e'5e-^'/' and r/_(r, z) = -2A^iCe<5e"^* - e'Je"^*/'^ 



arc respectively upper and lower solutions of the problem (2.8) provided e is sufficiently small. Hence, 
by the maximum principle we see that 

-2X-'CeSe-^* - e'Se'^*^^ < fi{r, z) < 2X-^CeSe-''' + e'Se'^'^^ (5.8) 

for ah < r < 1 and t > 0. 

Next, by differentiating the equation in the first line of (2.8) in r and using the boundary value 
conditions for rj, we see that fir satisfies 

ee^'^'\vr)t = {flr)rr + + ee^<'\^{7i,p, z)\{nr)r ' [e2<*)a(r; 7?, ^) + ^ - ee^<'^s{v,P,z)]rir 

-ee'^'''^*'^^{ri,p,z)[mrz{r\z) -rmrr{r\z) -mr{r\z)\ for < r < 1, t > 0, ^^'^^ 
r?^(0,i)=0, ?7rr(l,i) + [2 + ffe2^^*)g(j7,p,^)]r?^(l,t) = -£e2^(*)g(?7,p,2)m^(l;^) for t > 0, 
f)r{r,Q) = -n'oir) for < r < 1. 

If s is sufficiently small we have 

e2^Wa(r; r],z) + ^- £e2^Wg(7?,p, z) > e^(^*-'^°)co + 2 - CeSo > X. 
Moreover, using the assumptions on r],z and also using (5.7) and (5.8), we see that 

|e^^(*)ar(r; ry, z)fj + e^''(*)a^(r; t], z)fjr]r + ee'^^^*^g{r],p, z) [m„(r; z) - rmrr{r; z) - mr{r; z)] \ 

Since for e sufficiently small we also have 

2 + £e2^(*)g(77,p, z) > 1 and \ee^'^'^ s{v,P, z) [m,{l; z) - m,(l; z)] \ < C6ee-''\ 

again by using the maximum principle we conclude that if £ and s' are sufficiently small then 

-(1 + 2X-^)CeSe-^'* - e'Se'^*/^ < fir{r, ^) < (1 + 2A-^)C£^e-''* + Ce'5e-^^l' (5.10) 

for all < r < 1 and t>Q. Combining (5.8) and (5.10), and choosing eq further small when necessary 
so that X/e > for < £ < eq, we obtain (5.5). This proves Lemma 5.1. □ 

Lemma 5.2 Under the same conditions as in Lemma 5.1, there also holds the following estimate: 



L fl 



Ja 



e-'^^*-'\firr{r,s)\''r^drds < CeS^e-^"* (5.11) 
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for all t > 0. 

Proof. Prom the equation in the first three hnes in (5.9) we have 

2 

£{fir)t = e"^''(*V"^(r^7?rr)r + erg(77,p,2)7?rr - [a(r; z) + ^e"^''^*) - es{r],p,z)]fir 
-ar{r; r], z)fi - ar,{r; r], z)f]r]r - ss{v,P, z) [mrz{r; z) - rmrr{r; z) - mr{r; z)] . 
Multiply both sides of this equation with fjrr^ and next integrating with respect to r in [0, 1], we get 

=-6"^''^*' / \'nrrfr^dr + sg{r],p,z) / fjrrfjrr^dr - / a{r;r],z)\fjr\^r'^dr 
Jo ^ Jo Jo 

_2e-2^(*) / \fj^fdr + sg,{r],p,z) / \fir\'^r'^dr - / ar{r\r},z)fir)rr'^dr 
Jo Jo Jo 

- / ar,{r;rj,z)fir]rfirr'^dr - e^{ri,p,z) I [mrz{r;z) - rmrr{r\z) - mr{r;z)\firr'^dr. 
Jo Jo 



We have 



ss{v,P,z)J firrfjrr^dr =-^{r],p,z)\rir{l,t)\'^ - —^{r],p,z) |r?r|Vdr 
KCe^S^e-^^"^ (by (5.5) and (5.7)), 

/ a{r;ri,z)\fjr\'^r'^dr>co[ |r?r|Vdr, 
Jo Jo 

2e-2^(*) / \jir\^dr>0. 
Jo 

Moreover, by using (5.5), (5.7) and the boundedness of r],p, z (with bounds independent of special choice 
of these functions) we have 

s:siv,P,z) / IVrl'^r'^dr - / ar{r;r], z)fifirr'^dr - / an{r;r], z)fjr]rfirr'^dr 
Jo Jo Jo 

-s^{r],p,z) / [mrz{r;z)-rmrr{r;z)-mr{r;z)]firr^dr 
Jo 

It follows that 

~ f'\fjr\'r'dr<-Co f \f,rr?r^dr-c^ f \f,r\\^dr + Ce'5\-^^' . 
^ dt Jo Jo Jo 

Prom this inequality we easily deduce that if s is sufficiently small such that < £ < |^ , then 

\fj,\^r'dr+^ f [\-"-^^'-'^\firr{r,s)\^r^drds 
Jo ^ Jo Jo 

<e-^* \v'o{r)\^r^dr + Ce^S^e-^^"' < CJ^e-^"* + Cs^6^e-^^*. 
Jo 
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The estimate (5.11) now immediately follows. □ 



CoroUciry 5.3 Under the same conditions as in Lemma 5.1, for any < u < p. there exists 
corresponding constant C > such that the following estimate holds: 

f f e^^*\firr{r,t)\'^r^drdt<C5^. (5.12) 
Jo Jo 

Proof. Multiplying both sides of (5.11) with e^"^* and integrating with respect to t over [0,00), we 

get 

/ e^"'-^^*-')\nrr{r,s)\''rHrdsdt<Cs5^ e-2(''-)*dt = -f^. 

pOO pt pi 

/ / / e^'''-^^'-'^\f,rr{r,3)\''r^drdsdt 
Jo Jo Jo 

pOO pi pOO 

7 y (y e-^^'^-''^'dt)\vrr{r,s)\^r^drds 



Since 



e 

"2(^ 

we see that (5.12) follows. □ 



e f°° 

7 / / \rirr{r,s)\'^r'^drds, 

— ^^j Jo Jo 



Corollary 5.4 Under the same conditions as in Lemma 5.1, for any < u < jj. there exists 
corresponding constant C > such that the following estimate holds: 

e''''\fit{r,t)\^r''drdt<—. (5.13) 
JO ^ 

Proof. This follows from Corollary 5.3 combined with the equation (2.7) and the fact that max |??(r, t) | < 

0<r<l 

Cde-''* and max \fir{r, t)\ < CSe'^* for all t>0. □ 

0<r<l 

6 The proof of Theorem 1.1 

The proof of Theorem 1.1: Let ji* be the positive constant specified in Lemma 4.2, and arbitrarily 
choose a positive constant ji such that < /U < /k* and fix it. Let 6 and e' be positive constants which we 
shall specify later. We assume that the initial data {cQ{r),po{r),zo) of {c{r,t),p{r,t), z{t)) satisfy (1-21), 
(1.22) and the following conditions: 

max |co(r) - c*(r)| < e'5, sup |co(r) - c*'(r)| < e'5, (6.1) 

0<'-<l 0<r<l 

max \po{'r) —p*{r)\ < 5, sup r{l — r)\p'Q{r) — p* (r)| < 5, and 1 2:0 — -2*1 < (6.2) 
o<i-<i 0<r<l 

Let (y, II • II y) be the Banach space introduced in Section 2, and consider the following set 5 C F: 

S = {rj^Y: r,{r,0) = vo{r), \\v\\y<S}, 

where r]o{r) = co(r) — m(r, zq). Let r] G S and Vq = {po,Zo). It is clear that Vo £ Xq and, due to the 
conditions in (6.2), we have ||Vb||xo < ^- It follows by Lemma 4.3 that there exists positive constant 
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5o depending only on the choice of /x such that ii < S < Sq then the equation (4.2) subject to the 
initial value condition V{0) = Vq has a unique solution V G C([0, oo),Xo) n C^([0, oo),X) satisfying the 
estimates in (4.18). Let vlt) = {ip{-, t)X{t)) and set p{r, t) = p*{r) + ip{r, t), z{t) = z* + ((t). Prom the 
deduction in Section 4 we see that {p,z) is a solution of the system (2.7). Clearly, {p,z) satisfies (1.26) 
and (1.27), so that it also satisfies (5.2) and (5.3) when we replace 5 in those conditions with CS. Now 
we consider the problem (2.8). It is clear that the condition (5.1) is satisfied and, by (6.1) and (6.2), the 
condition (5.4) is also satisfied. It follows by Lemma 5.1 that there exist constants eo, Sq > independent 
of S such that if < £ < £o and < e' < e'q then the solution fy = ri{r,t) of (2.8) satisfies the estimate 
(5.5), so that 

\fj{r,t)\ + \fir{r,t)\ < C{s + e')6e-^'' < ^e""* 

provided C(eo +eo) < 1- This implies that fj G S ii e and e' are sufficiently small. Let ^ : S ^ S he the 
map ^{r]) = t). In what follows we prove that ^ has a fixed point. 

We first prove that ^ is continuous with respect to the metric 

rfi(??i,??2) = supe'^*f / \'ni{r,t)-'n2{r,t)\^r^drY . 
Given r/i, 772 € 5, we let f\i = ^{r]i), i = 1,2. Then fji is the solution of the following problem: 

efjit = e~'^''*^*\~'^{r'^fjir)r + £rg{r]i,pi, Zi)fjir - a{r;r]i, Zi)fji 

-eg{r]i,pi, Zi) [rrizir; Zi) - rmr{r; Zi)] for < r < 1, t > 0, 
77,^(0, t) =0, J7,(l,t) =0 for t>0, 
f]^{r, 0) = ?7o(r) for < r < 1, 

where {pi,Zi) is the solution of the problem (2.7) for rj = rji, i = 1,2. It follows that 

.1 .1 

Jo 

-ysCm^Pi^^^i) J |^i-^2|Vdr + e[g(?7i,pi,2;i)-g(?72,P2,^2)] J [vi - V2]mrr^dr 



=_e-2^i(*) 



a{r;r]i,zi)\iii - fj2fr'^dr - / [a{r;r]i, zi) - a{r;r]2, Z2)][fii - mjmr'^dr 
1 

/ {&{m,Pi,zi)[m;,{r;zi) -rmr{r;zi)] - g{'n2,P2, Z2)[m;,{r; Z2) - rmr{r; Z2)]} 
Jo 

X [vi — fiiV^dr. 

Since a(r; 771,21) > cq > and rji, pt, Zi {i = 1,2) are bounded functions with bounds independent of 
special choice of these functions, by using (5.5), (5.7) and some standard arguments we see that if 6 is 
sufficiently small then 



-2zi(t) 



\flir -mA r dr - -Co 



\m - 772! r dr 



+CS^e-^^*[max \r]i{r,t) - r]2{r,t)\^ + max \pi{r,t) - P2{r,t)\^ + \zi{t) - Z2{t)\^] 

0<r<l 0<r<l 

+Ce\m^x \m{r,t) - V2{r,t)\^ + max \p,{r,t) - P2{r,t)\^ + \z^{t) - ^2(i)l']- 

0<r<l 0<r<l 
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Neglecting the first term on the right-hand side of the above inequality, assuming that e is sufficiently 
small so that Ce < ^cq, and using the estimate (4.26), we get 

~ /'|7?i-7?2|Vdr<-^ f \f,^-f,^\\^dr + C{5^ + e^)e-^>^' sup e''^'\m{r,t) ~ m{r,t)\\ 
I at Jq z Jq o<r<i 

t>0 

It follows that by assuming that e is further small such that cq — 2/i£ > ^cq when necessary, we have 

[' \m{r,t)-m{r,t)\^r^dr<C{S^+s^)e-^>''( sup e''*|r?i(r,i) - ?72(r,i)l)'- (6.3) 

Jo ^ 0<r<l / 

t>0 

By applying the three-dimensional interpolation inequality 

||w||oo <C||Vw||i||M||| for u G H^{Bi{0)), Vu G L°°{Bi{0)) 
to the case u{x) = u{\x\), we get the following inequality: 

3 />! 1 

sup \u{r)\<c( sup \u'{r)\y( \u{r)fr^dr)\ (6.4) 

0<r<l ^0<r<l ^Jo ' 

It follows that for any , S we have 

sup e''*|r?i(r,i) - r]i{r,t)\ < CSi sup fe"'^* / |r/i(r,i) - r]2{r,t)\^r^drY . (6.5) 

o<T-<i (>0 ^ Jo ' 

t>o — 

Substituting (6.5) into (6.3) we get 

Hence, the map ^ : 5 — )■ 5 is continuous with respect to the metric d\. 
Let 

5'o =1'? G 5 : r]{r,t) is twice weakly differentiable in r and weakly differentiable in t, 

/ / e'**|r?^^(r,t)|Vrfrdt<C<52 and / / d'%t{r,t)'^r'^drdt <—\, 
Jo Jo Jo Jo ^ ^ 

where C is the constant appearing in (5.12) and (5.13) for the case v = /x/2. By Lemma 5.3 and Lemma 
5.4 we see that ^(5) C 5o. 

Let S be the closure of S with respect to the metric di , and let : S* S' be the unique continuous 
extension of the map ^ onto S. Then ^{S) C So, where So is the closure of with respect to the 
metric di . By using a standard *-weak compactness argument we easily see that 

-So c{r/ e C{[0,oo),L'^((Q,^),r^dr))nL^((0,oo),H^{{0,l),r^dr),e^''dt) 

nffi((0,oo),L2((0, l),r2dr),e'^*dt) : r]{r,0) = rio{r) for a. e. r G (0, 1), 
r]{l,t) = for a. e. t > 0, esssup |77(r,t)| < 6, esssup e'^ |ryr-(r,t)| < S, 

0<r<l 0<r<l 
f>0 t>0 

' / e^*|r?^,(r,t)|V2cirdi < CJ^ and / / e^'\T]t{r,t)\^r' drdt < [. 

Jo Jo Jo ^ ^ 
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Let 



(■1 1 

F= {r?e C([0,oo),L2((o,l),r2dr)) : ||77||i. =^ supe"* ( j \r]{r,t)\''r''dry < oo}. 



It is clear that (V, \\ ■ ||-p) is a Banach space and 5 is a closed convex subset of this space. For any 
T > 0, wc let Yt, Yt, St, St, Sqt and Sqt be respectively the restrictions of Y, Y, S, S, So and S^o on 
[0,1] X [0,T], and correspondingly define ,^t, -^t to be respectively the "restrictions" of ^ on St 
and St, i-c, for any G St we choose an 77 e 5 such that r?|[o,T] = C and define ^t{0 — ■^{i)\[o.t] , and 
similarly for By uniqueness of solutions of the problems (2.7) and (2.8) in any interval [0,T] wc see 
that these definitions make sense. Moreover, it is clear that ^t is a continuous self-mapping in St, ^t 
is a continuous self-mapping in St, .'^t{St) Q Sqt and ^t{St) Q Sqt- Note that, in particular, 

Yt = {v& C{[0,T],L\{0,l)ydr)) : My^ sup e^*( |r?(r, t)| V^dr) ' < ooj, 

0<t<T ^Jo ' ' 

SotC {v e C{[0, T], L\{0, 1), r^dr)) n L^{0, T), H^O, 1), r^dr)) n H'{{0, T), ^^((o, l)ydr)) : 
r]{r,0) = r]o{r) for a. e. r e (0, 1), r]{l,t) = for a. e.t€ (0,T), sup e'^*\r){r,t)\ < S, 



0<r<l 

0<t<T 

fT /.I pT 1-1 



sup e'^* I r/^(r,t) I < (5, / / e''* | r/^^ (r, t)p r^drdt < CJ^ and / / e>'^\j]t{r,t)\^r'^drdt < }. 

0<'-<i Jo io Jo io ^ 

Clearly, St is a closed convex subset of Yt- Moreover, it is also easy to see that Sqt is a compact 
subset of Yt- Indeed, let {?7n}^i be a sequence in Sqt- By compactness of the embedding _ff^((0, 1) x 
(0, T), r^drdt) ^ L^qq, 1] x [0, T],r'^drdt), *-weak compactness of bounded sets in L°°([0, 1] x [0, T]), and 
weak compactness of bounded sets in i^([0, 1] x [0, T]), we see that there exists a subsequence of {?7rt}^i 
which we assume, for simplicity of the notation, to be {?7n}^i itself, and r] G Sqt such that as n — )• 00, 

r]„^r] strongly in L^([0, 1] x [0, T], r^drdt), 

Vrn Vr *-weakly in L°°([0, 1] x [0, T]), 
Vrrn " > Vrr wcakly in L^([0,1] X [0,T],r'^drdt), 
Vtn Vt weakly in L'^{[0, 1] x [0, T],r'^drdt). 

Integrating the equation 
dt , 



■ / \r]n{r,t)-r]{r,t)\ydr = 2 [ [r]n{r,t) - r]{r,t)][r]nt{r,t) - vM^dr, 
Jo Jo 

and using the Cauchy-Schwartz inequality and the inequalities / / e^*\rjnt{r,t)\'^r'^drdt < (n = 

Jo Jo ^ 

l,2,---)and / / e^*'\r]t{r,t)\'^r^drdt < , we get 

Jo Jo ^ 

sup / \T]n{r,t)-r,{r,t)\''r^dr = ^( / \7jnir,t) - f]ir,t)\^ydrdt)\ n = l,2,---. 
0<t<TJo V£ ^ Jo Jo ^ 



It follows that lim ||r/„ — tjWy^ = 0. This proves the desired assertion. 



T 

V ■ 



It follows by the Schauder fixed point theorem that ^t has a fixed point in St, which we denote as 
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We now prove that ry'" G Sqt and it is a fixed point of We first prove that rj^ G C'([0, 1] x [0, T]). 
This assertion follows from the fact that Sqt C C([0, 1] x [0,T]), which is proved as follows: First, for 
any p > 1 by integrating the equation 

\'n{r,t)\^r^dr =p j \r]{r,t)\P~'^sgnr]{r,t)r]t{r,t)r'^dr 

and using a similar argument as before we see that for any rj G Sqt such that 0) = 0, we have 

sup r \r]{r,t)\Pr'^dr = 9^lJ-f \r]{r,t)\^r^drdtY . 

<t<TJo v£ ^Jo Jo ' 



0<t<T Jo 

Next, by integrating the equation 



^\v{r,t)\^ =2n{r,t)r]r{r,t) 



and using the Holder inequality for p > 3 we see that for any r] G Sqt we have 



sup |?7(r,i)| <2 sup / \ri{r,t)\dr ■ sup Irjrifjt)] < C6 sup / \r]{r,t)\rp ■ r pdr 

o<T<i 0<t<TJo o<r<i 0<t<TJo 

0<t<T ~ ~ 0<t<T ~ ~ 

' \r]{r,t)\Pr'^dr)'' 
^ 



Combining the above two inequalities we get 



sup \r,{r,t)\ < ( / / \j]ir,t)\''r''drdt)"' . (6.6) 



0<'-<i ^ ^Jo Jo 

0<t<T 



Using this inequality to ry„ — r/ for any r] G Sqt and a corresponding sequence r/„ e Sot {n = 1, 2, • • • ) 
such that lim / / |r7„ — r^Pr^drrft = 0, we see that r7„ converges to uniformly in [0, 1] x [0, T], so that 

T] e C([0, 1] X [0,T]). This proves that Sqt Q C([0, 1] x [0,T]) and, consequently, r;^ G C([0, 1] x [0,T]). 

It follows that the problem (2.7) with rj replaced by ry^ has a unique classical solution for < t <T, 
which we denotes as (p'^,z'^). Since rf ^p^ G C([0, 1] x [0,r]) and G C[0,T], by using the standard 
L^-theory for parabolic equations we see that when (?7,p, z) is replaced by {rf ,p^ , z^), the problem (2.8) 
has a unique strong solution in [0, 1] x [0,T] which we denote as ff . Clearly, G Wp'-^((0, 1) x (0,T)) 
for any 1 < p < oo. Using a limit argument we can easily show that fj^ = rf^ . Indeed, by choosing a 

sequence r/„ G Sqt {n = 1,2,- ■ ■) such that lim / / |ry„ — ry^pr^rfrdt = (which implies, by (6.6), that 

Jo Jo 

also lim sup |ry„ — 'rf\ =0) and defining r)„ = ^riVn), we see that ^„ — )• fj^ uniformly in [0, 1] x [0, T]. 

rn-oo o<r<l 

0<t<T 

Since also ^„ = ^riVn) ^rif}^) = f]^ strongly in Yr, by uniqueness of the limit we obtain the desired 
assertion. It follows that (rj'^ ,p'^ , z'^) is a solution of the problem (2.6) in [0, 1] x [0,T]. Using this fact 
and the bootstrap argument we immediately obtain the assertion that tj^ G C([0, 1] x [0,T]). Hence 
T]^ G Sot- The assertion that r]^ is a fixed point of ^t now follows from the fact that r)-^ = r]'^. 

Since {rf" ,p'^ , z'^) is a classical solution of the problem (2.6) in the time interval [0, T], by imiqueness 
of the solution of this problem, we conclude that for any < Ti < r2 we have rj'^^ (r, t) — rj^'^ (r, t) for 
(r, f) G [0,1] X [0,ri]. It follows that the following definition of the function rj in [0,1] x [0, oo) makes 
sense: 

7?(r, t) = ri^{r, t) for any (r, t) G [0, 1] x [0, T] and T > 0. 
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Moreover, by letting T — > oo in the relations ry^ e Sqt and ^riv'^) = , we see that r] G So C S and it 
is a fixed point of the map : S ^ S. 

Having proved that ^ has a fixed point in S, the desired assertions in Theorem 1.1 then immediately 
follow from Lemma 4.3 and Lemma 5.1. This completes the proof of Theorem 1.1. □ 
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